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ABSTRACT 

The  subject  of  this  paper  is  the  buckling  behavior  of  a  rectangular  plate,  with  parallel  thin- 
walled  stiffeners  attached  to  one  side,  subjected  to  a  combination  of  axial  compression,  latero.l 
pressure,  and  bending  moment.  The  plate  is  modeled  by  the  Von  Kdrmdn  plate  equations  and 
the  stiffeners  by  a  nonlinear  beam  theory  recently  derived.  An  analytical  solution  is  obtained 
for  the  buckling  load  corresponding  to  a  torsional  tripping  mode  of  the  stiffeners.  The  effects 
of  various  boundary  conditions,  imperfections,  and  residual  st\t^s  are  included. 
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1.  INTRODUCTION 

Stiflfened  plates  are  a  basic  structural  component  of  ships  and  submarines.  These  struc¬ 
tures  are  designed  with  generous  safety  margins  against  overall  collapse  triggered  by  buckling. 
The  object  of  analjrtical  work  is  to  determine  design  criteria  to  inhibit  buckling  at  any  stress 
less  than  yield.  Recently  [see  Danielson  et  a/  (1993)],  we  have  developed  an  analytical  formula 
for  the  buckling  load  of  a  stiffened  plate  subjected  to  a  combination  of  axial  compression 
and  lateral  pressure.  The  object  of  the  present  paper  is  to  improve  and  extend  our  previous 
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analysis.  A  review  of  the  literature,  given  in  our  earlier  work,  will  not  be  repeated  here  but 
details  of  the  analysis,  which  supersedes  our  earlier  work,  will  be  recorded  here. 

We  Srst  consider  a  plate  which  is  initially  rectangular  in  shape  and  has  several  parallel 
I-stiffeners  spaced  a  distance  b  apart.  The  structure  is  subjected  to  a  combination  of  uniform 
axial  compressive  stress  a  (force  per  unit  area  of  a  side),  uniform  lateral  pressure  p  (force  per 
unit  lateral  area  of  the  plate),  and  uniform  bending  moment  M  (moment  per  unit  length  of 
an  edge).  We  suppose  that  at  low  values  of  cr,p,  and  M  the  plate  and  stiffeners  simply  bend 
and  compress  symmetrically.  Our  object  is  to  find  the  critical  load  at  which  the  stiffened 
plate  may  buckle  into  an  alternate  mode  (see  Fig.  1). 

Our  present  analysis  is  based  on  the  following  simplifying  assumptions: 

(i)  Each  plate-stiffener  unit  of  width  b  undergoes  an  identical  deformation. 

(ii)  The  plate  obeys  the  honlinear  Von  Karman  plate  equations  [see  Timoshenko  and  Gere 

(1961)].  The  stiffeners  obey  the  nonlinear  beam  equations  derived  by  Danielson  and 
Hodges  (1988). 

(iii)  The  plate  and  stiffener  material  is  elastic,  linear,  and  isotropic. 

(iv)  Every  particle  on  the  bottom  surface  of  a  beam  undergoes  the  same  displacement  as 
the  corresponding  particle  on  the  top  surface  of  the  plate,  and  every  line  of  particles 
in  the  beams  normal  to  the  plate  surface  remains  normal  to  the  deformed  plate  at  its 
surface.  In  other  words,  the  bases  of  the  stiffeners  are  clamped  to  the  plate. 

(v)  The  prebuckling  displacements  are  less  than  the  maximum  thickness  of  the  structure 

and  independent  of  the  transverse  coordinate. 

(vi)  The  incremental  buckling  extensional  strains  at  the  midsurface  of  the  plate  are  negli¬ 
gible. 
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(vii)  The  incremental  buckling  displacements  may  be  approximated  by  the  fundamental 
harmonic  in  their  Fourier  expansions. 

(viii)  The  plate  and  beams  are  so  thin  that  their  thicknesses  are  negligible  compared  to 
their  width,  height,  length,  and  the  wavelength  of  deformation.  A  stiffener  is  so  slender 
that  its  width  and  height  are  negligible  compared  to  its  length  and  the  wavelength  of 
deformation. 

2.  GENERAL  POTENTIAL  ENERGY  FUNCTIONAL 

It  follows  from  assumption  (i)  that  we  need  only  analyze  a  single  plate  unit  containing 
a  single  stiffener.  From  a£;sumptions  (ii)-(iii),  the  potential  energy  of  the  plate  plus  beam  is 
given  by: 

/4  f(x  •  P* f 

Jo  ^  +  (1  -  i/)ey 

+£>  ^  +  (1  -  ^')u^?2|  +  —  P^  +  Aftynj  dxi  dx^ 

^  J  I  Jb  +  2G'|?3)  ^^2  dx3  +  (tA{ui  -  cwn)dxi 

(1) 

Here  (xi,  xj,  X3)  are  Cartesian  coordinates  measured  from  the  midpoint  of  a  side  of  the  plate. 
The  plate  unit  has  length  a,  width  6,  and  thickness  f,  while  the  beam  has  cross-sectional 
area  A  and  centroidal  height  c.  The  elastic  constants  are  defined  by 

n  a--I— 

12(1 -./«)’  2(1  +  ./) 

where  E  is  the  Young  modulus  and  1/  is  the  Poisson  ratio.  The  displacements  of  the  plate 
midsurface  in  the  xi,X2,  and  X3  directions  are  denoted  by  u(xi,  xj),  v(xi,X2),  and  iu(Xi,  X2), 
respectively.  Subscripts  on  u,i;,  or  w  denote  partial  differentiation  with  respect  to  the 
coordinates  Xj  or  xj;  e.g.,  W12  =  The  extensional  strains  at  the  midsurface  of  the 

plate  are  given  by 
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1  2  1  1 
Cll  -  Ui  +  -W, ,  Ci2  =  -(U2  +  +  U?iU;2),  622  =  1’2  +  ^“^2 

The  strains  in  the  beam  are  denoted  by  7ii(xi,  ij,  X3),7i2(®i,  X2,  X3),  and  713(2:1,  X2,  2:3) 
and  are  related  to  the  displacements  by  equations  (9,  10)  of  Danielson  and  Hodges  (1988), 
which  upon  invoking  assumption  (iv)  are  transformed  into: 

711  =  ■£'11  +  (■£'12  +  £12)^3  “  j£'13<^2  +  2*^2  "h  2^3 

712  =  /?12  —  7l3  =  ^13  +  ■^E\\(f>2 

En  =  Cll  —  Xsuiii  +  Aujii2,  E\2  =  612  —  2^3^12  +  2^2^12 

^  1  _  1  ,  _  1 , _ 

■£'13  —  2^21^12  +  2'^3li’l2  + 

,  1  _  1  ,  _  1  , _ 

<t>2  —  “2^2^12  +  2''3t^l2  “ 

<i>3  =  —  2^3U'12  —  2'^2^12 ' 

Here  A(x2i  2:3)  is  the  Saint- Venant  warping  function  for  the  beam  cross  section;  subscripts 
on  A  denote  partial  differentiation  with  respect  to  X2  and  X3.  Bars  over  a  symbol  denote  its 
value  at  the  beam  axis;  e.g.,  the  axial  displacement  of  the  beam  centroid  is  u(xi)  —  dl;i(xi). 

Substituting  these  relations  into  (1)  and  neglecting  higher  than  cubic  terms  in  the  dis¬ 
placements  (these  are  not  needed  in  our  subsequent  analysis),  we  obtain  a  lengthy  expression 
for  the  piotential  energy  which  forms  the  basis  for  our  subsequent  analysis.  Among  all  the 
functions  satisfying  the  geometric  or  natural  boundary  conditions  the  one  which  causes  the 
potential  energy  to  be  a  minimum  is  the  equilibrium  state.  We  suppose  that  the  outer  edges 
of  the  plate  arc  free  to  displace  in  the  horizontal  plane  but  restrained  in  the  vertical  direction, 
so  the  geometric  boundary  conditions  are 

u;(0,X2)  =  u;(a,X2)  =  0  (2) 

Note  that  the  case  of  simply  supported  edges  is  obtained  by  setting  A/  =  0. 
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3.  PREBUCKLING  SOLUTION 


The  prebuckling  equilibrium  state  is  denoted  by  (u,w).  It  follows  from  assumptions 
(ii)-(v)  that  the  potential  energy  in  the  prebuckling  state  is 

Jo 

Here  I22  is  the  moment  of  inertia  of  tb?  beam  section  about  the  xj-axis: 


2  2 

+cr(i4  +  ib)  ui  —crAc  wn  —pb  w  +Mb  u^uj  dxi 


h7  =  /  3^3  <^^2  dx:^ 


The  prebuckling  displacements  are  determined  by  the  variational  equation 


8P  =  0 


Taking  the  variation  of  (3),  and  integrating  by  parts  with  respect  to  xj,  we  obtain  the 


differential  equations 


{bD  +  EI22)  njiiii  —EAc  Uju  —pb  —  0 
—EAc  u;ii  +E{A  +  tb)  Uj  ■\-ff{A  +  ft)  =  0 


and  boundary  conditions 


at  X  =  0  and  x  =  a  :  {bD  +  EIzi)  ^11  —EAc  Uj  —a Ac  +  Mb  —  0 


The  solution  to  the  linear  boundary  value  problem  (2),  (4),  (5)  is 


o  (J  Ac  u;ii 

“  E^  A -\-tb 


e 

W  = 


px,(a  -  xi)(a^  +  gxi  -  x\) 


Mxi{a  —  xj) 


Note  that  (6)- (7)  reduce  in  the  case  /jj  =  ^4  =  0  to  the  well-known  exact  solution  for  an 
isolated  wide  plate,  and  in  the  case  f  =  0  to  the  well-known  exact  solution  for  an  isolated 
beam. 
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4.  BUCKLING  SOLUTION 


According  lo  llie  energy  criterion  of  elastic  stability,  the  prebuckllng  equilibrium  slate  Is 
stable  if  and  only  if  the  energy  functional  which  represents  the  increase  of  the  total  potential 
energy  in  a  displacement  field  lo  some  slightly  adjacent  slate  (u;  +u;)  is  non-negative: 

P[^,+u,]-P[u;]>0.  (8) 

Since  the  prebuckling  state  is  an  equilibrium  state,  the  terms  in  (8)  which  are  linear  in  the 
incremental  displacement  w  must  vanish.  It  follow's  that  the  terms  Q[tu]  in  (8)  which  are 
quadratic  in  the  incremental  displacement  must  be  non-negative: 

Q(uj]  >  0 . 

The  critical  case  of  neutral  equilibrium  occurs  when  there  exists  a  buckling  mode  Wcr  satis¬ 
fying 

Qkcr]  =  0  (9) 

C?[it;  5^  ujcr]  >  0.  (10) 


The  eigenvalues  (Ter, Per,  and  Mcr  which  render  (9)  zero  are  the  critical  buckling  loads. 

From  the  first  integral  in  (1)  and  assumption  (vi),  the  quadratic  functional  for  the  plate 


+  l/WnW22  +  —  4- {1  -  t/)u!l2  - 


dx\dx2.  (11) 


From  the  remaining  integrals  in  (1),  the  quadratic  functional  for  the  beam  is 


L  2  h 


*  -2-  ^2 
«  It  _  _  Aft* 


tyji  lujjdxi  + 


{GJ  -  al)  r. 


f  njJjdxi . 

Jo 


Here  /  is  the  polar  moment  of  inertia  about  the  xi-axis,  J  is  the  Saint- Venant  torsion 
constant,  aind  Hi,  H2,  H3  are  constants  defined  by  the  following  integrals  over  the  beam 
cross  section: 


^  I  /bom  (^2  ^^3 

J  =  I  L,t„  [(2:2  +  Aa)^  +  (13  —  A2)^]<fx2  dx3 
=  /  L,„,  X^dxidxz 

H2  =  I  L^„[x3{xl  +  xl- Xl-  Xl)  +  2\{X2  + X3)]dx2dX3 

J  •/••cliea 

H3  =  7  /  /b„„  a:3[(x2  -  Aa)^  +  (xa  +  A2)’  +  2(x2  +  X3  -  A^  -  A3)]<fx2  dx3 

4  J  t/MClion 

The  total  quadratic  functional  for  the  plate  plus  the  beam  is  the  sum  of  (10)  and  (11). 

Next,  we  calculate  the  cross  section  properties  for  a  beam  composed  of  a  thin  web  and  a 

thin  bottom  and  top  flange.  The  web  has  thickness  and  height  hru\  the  bottom  flange  has 

thickness  tj  and  width  /ij,;  the  top  flange  has  thickness  if  and  width  hj.  The  Saint- Venant 

warping  function  for  this  thin-walled  cross  section  is: 

{X2{2h^  +  tf  -  X3)  flange 
X2X3  web 

-X2X3  bottom  flange 

Using  approximation  (viii)  we  obtain: 


=  h^b  +  t  +  </  h 

__  ,  J  l2  L 

=  K 

t;*.  iiK 
3  ^  3  ^  3 
tjhlh} 

12 

tjhl 

3  6 

=  -y -b  t/ h/ +  . 
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In  accordance  with  approximation  (vii)  and  the  boundary  conditions  (2)  and  (5),  we 
represent  the  incremental  buckling  displacement  by  the  following  shape  (an  arbitrary  multi¬ 
plicative  constant  has  been  set  equal  to  1): 

(13) 


.  VlTTXi  .  irX2 

10  =  sin - sm— m  =  l,2,3 

a  0 


Substitution  of  this  buckling  mode  into  (11)  plus  (12)  and  application  of  the  inequality  (10) 
leads  finally  to: 


CTcr  < 


F 


„  E,,  , 


vjLy  .  CL  .9  T^m^EH,  „ , 


bD  .mb 
2  '  a 


r 

^  +  57 


(14) 


Here  m  is  taken  to  be  the  integer  which  gives  the  lowest  value  of  <Tcr  in  (14).  Note  that 
(14)  reduces  in  the  case  h  —  iu,  =  ij  =  0  to  the  well-known  exact  solution  for  an  isolated 
plate,  in  the  case  per  =  M„  =  i  =  0  lo  our  previous  exact  solution  for  an  isolated  beam  [see 
formula  (28)  of  Danielson  et  al  (1990)],  and  in  the  case  pcr  =  ^ct  =  0  to  formula  (87)  (with 
=  (<^e)^e  and  -k  F  w  Hi)  of  Adamchak  (1979). 


5.  CLAMPED  EDGES 


In  this  section  we  consider  the  clamped  case  when  the  bending  moment  M  is  not  pre¬ 
scribed,  but  the  rotation  at  the  edges  is  completely  restrained.  Then  the  additional  geometric 
boundary  conditions  are 

iui(0,X2)  =  uj,(a,X2)  =  0  (15) 


Note  that  in  this  case  each  material  peu'ticle  on  an  edge  cross-section  is  totally  restrained 
frcm  any  motion  at  the  buckling  point. 

The  prebuckling  solution  to  the  linear  boundary  value  problem  (2),  (4),  (15)  is  (6)  and 

px]{a  -  xif 


o 

W= 


F 

241D  +  -(/„-^^)1 


(16) 
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Note  that  (6)  and  (16)  reduce  in  the  case  7^2  =  /I  =  0  tc  the  well-known  exact  solution  for 
an  isolated  wide  plate,  and  in  the  case  <  =  0  to  the  well-known  exact  solution  for  an  isolated 
beam. 

In  accordance  with  approximation  (vii)  and  the  boundary  conditions  (2)  and  (15),  we 
now  represent  the  incremental  buckling  displacement  by  the  following  shape: 


mTTXi  .  TTIj  „  .  . 

10  =  (1  —  cos - )sm— j— ,  m  =  2,4,6,---  (17) 

a  b 

Substitution  of  this  buckling  mode  into  (11)  plus  (12)  and  application  of  the  inequality  (10) 
leads  finally  to 

I  Ac 

^  bD..mb  ,  a  .  20^  TrWEHj  .  ^  , 


< 


4ir!m=lZ)  +  -  -f— j)] 


A  tb 


(18) 


Note  that  (18)  reduces  in  the  case  per  =  i  =  0  to  our  previous  exact  solution  for  an  isolated 
beam  [see  formula  (28)  of  Danielson  et  al  (1990)]. 


6.  IMPERFECTIONS 

In  this  section  we  suppose  that  the  structure  has  an  initial  normal  deflection  in  the  shape 
of  the  prebuckling  normal  displacement  caused  by  the  pressure  p.  Specifically,  when  M  is 
prescribed  on  the  edges  the  initial  normal  deflection  is 

16xi(a  -  ii)(a’  +  ax,  -  xj) 

- 5^5 - 

while  for  damped  edges  the  initial  normal  deflection  is 

We  also  assume  that  the  amplitude  W  of  the  initial  displacement  is  less  than  the  maxi¬ 
mum  thickness  of  the  structure.  Then  the  prebuckling  displacements  are  still  given  by  our 
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previous  solutions  (6)-(7)  or  (16),  and  the  quadratic  functional  Qpute  >s  still  given  by  (11). 
The  only  effect  of  this  prebuckling  deflection  is  to  create  a  new  term  in  Qbe»m  which  is  the 
same  as  the  naiddle  integral  in  (12)  with  the  prebuckling  normal  displacement  replaced  by 
the  initial  displacement,  and  we  can  use  our  previous  calculations  to  evaluate  this  integral. 
We  thereby  find  that  the  effect  of  this  imperfection  is  to  add  an  additional  term  to  our 


=  <r„llV  =  01  + - -  (20) 

Note  that  for  an  asymmetrical  structure  torsional  deformation  of  the  stiffeners  may  ini¬ 
tiate  upon  application  of  the  slightest  load,  so  bifurcation  may  not  be  able  to  be  used  as  the 
buckling  criteria  [see  Ostapenko  and  Yoo  (1988)]. 


7.  RESIDUAL  STRESS 

The  simplest  way  to  account  for  residual  stress  is  to  assume  that  the  plate  is  subject  to 
a  unifOTm  compressive  residual  stress  S,  while  the  beam  is  subject  to  a  counterbalancing 
distributkm  of  residual  stress  <T,(x2,i3)  [see  Hughes  (1983)].  The  only  effect  of  this  residual 
ftrcas  B  to  create  a  new  term  in  Qp,..,  which  is  the  same  as  the  last  integral  in  (11)  with 
o  replaced  5.  Note  that  the  analogous  term  in  Qbeam  is  zero  because  uj(xi )  =  0  for  the 
assoined  mode  shapes  (13)  and  (17): 

~  I J  L  ~  ^ 
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We  thereby  find  that  the  effect  of  this  residual  stress  is  to  add  an  additional  term  to  our 
formula  for  <Tcr' 

S 


Oct  -  =  0]  - 


1  + 


2^2/ 


(21) 


CONCLUSIONS 


Simple  analytical  formulas  that  include  the  effects  of  combined  loading,  various  boundary 

conditions,  imperfections,  and  residual  stress  do  not  appear  to  exist  in  the  literature.  As 

an  example  of  the  numerical  predictions  of  our  formulas,  let  us  assume  the  following  typical 

parameter  values  [taken  from  Smith  (1975)]: 

E  =  30,000  ksi  h  =  0 

1/  =  .3  /ifc  =  0 

a  =  48  in  =  .28  in 

6  =  24  in  =  5.5  in 

t  =  .31  in  if  —  .56  in 

hf  —  3.1  in 

For  this  example,  the  formulas  (14),  (19),  (21)  and  (18),  (20),  (21)  reduce  to 
M  case  (m  =  1): 

cr„  =  49  +  60pcr  +  ■45Mcr  +  55W  -  .765  (22) 

Clamped  case  (m  =  2): 

<rcr  =  105  +  6.6pcr  +  30W  -  .765  (23) 

The  collapse  loads  of  Ccr  =  27.8  psi  (per  =  0)  and  <7cr  =  27.1  psi  (pcr  =  .015  ksi)  measured 
by  Smith  (1975)  on  ship  grillages  may  be  accounted  for  by  choosing  appropriate  values  of 
A/er,  W,  and  5  in  (22)  or  (23). 

We  have  made  an  attempt  to  verify  the  accuracy  of  some  of  the  approximations  upon 
whidi  our  analysis  is  based.  For  instance,  in  the  M  case,  we  included  the  effect  of  beam 
cross-sectional  deformation  by  allowing  the  web  to  undergo  a  lateral  buckling  displacement 

V  =  sin  -I-  C2X3) , 
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where  Oi  and  Cj  are  constants  determined  by  minimizing  the  critical  axial  stress  o’er-  For 
the  parameter  values  listed  above,  this  effect  turned  out  to  be  negligible.  It  is  possible  to 
invent  unusual  cases  in  which  the  cross*sectional  deformation  is  of  importance,  but  for  most 
practical  dimensions  the  assumption  of  a  rigid  cross-section  seems  ok. 

For  another  instance,  in  the  clamped  case,  we  added  together  the  m  =  2  and  m  =  4 
normal  buckling  displacements  (17): 

.  2irxi  4nxi  ,  .  7rx2 

ty  =  [1  —  cos - f-  0/3(1  —  cos - )]  sm  — j—  , 

a  a  0 

where  C3  is  a  constant  determined  by  minimizing  the  critical  axial  stress  cTct-  For  the 
parameter  values  listed  above,  this  effect  also  turned  out  to  be  negligible.  For  most  practical 
dimensions  the  assumption  of  a  simple  buckling  mode  seems  ok. 

At  any  rate,  if  we  include  the  additional  displacement  functions  needed  for  the  above 
effects,  it  doesn’t  seem  possible  to  obtain  simple  formulas  for  the  unknown  coefficients  or 
buckling  loads. 
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Figure:  STIFFENED  PLATE,  LOADING  APPLIED,  BUCKLING  MODE 
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